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LONG-TERM  GOALS 

The  long-term  goals  of  “Radiance  in  a  Dynamic  Ocean  (RaDyO):  Radiance  and  Visibility  as  Affected 
by  Inherent  Optical  Properties”  (award  number  N0001409C0443)  and  “Imaging  System  Performance 
and  Visibility  as  Affected  by  the  Physical  Environment”  were  to  perfonn  detailed  investigations  of 
inherent  optical  properties  (IOPs)  and  the  effects  of  IOP  variability  on  underwater  radiance  and 
visibility  and  to  investigate  the  relationships  between  physical  forcing  and  optical  properties,  including 
the  imaging  performance  parameters.  These  efforts  were  aimed  at  supporting  the  Radiance  in  a 
Dynamic  Ocean  (RaDyO)  program,  which  was  devoted  to  the  topic  of  light  propagation  and  imaging 
across  the  air-sea  interface  and  within  the  surface  boundary  layer  of  natural  water  bodies. 

OBJECTIVES 

The  objectives  of  these  RaDyO  projects  were  to: 

(1)  Investigate  the  differences  between  optical  properties  and  optical  variability  measured  off  the 
R/P  FLIP  and  off  the  R/V  Kilo  Moana  (KM)  during  the  RaDyO  Santa  Barbara  Channel  (SBC) 
experiment  and  the  effects  of  optical  variability  on  the  modeling  of  the  underwater  radiance 
distribution  and  visibility. 

(2)  Detennine  the  relationships  between  physical  processes  and  the  optical  properties  that  affect 
underwater  visibility  and  the  imaging  perfonnance  parameters. 

(3)  Investigate  the  effects  of  vertical  variability  (e.g.,  vertical  layering)  of  optical  properties  on  the 
image  performance  parameters,  i.e.  modulation  transfer  function  (MTF)  and  point  spread 
function  (PSF). 

APPROACH 

Objective  (1) 

Concurrent  profiles  of  optical  properties  collected  off  the  FLIP  and  the  KM  during  the  SBC 
experiment  were  analyzed  for  optical  variability.  The  datasets  from  the  two  platforms  were  merged 
with  depth  and  density  and  compared.  Percent  differences  between  KM  and  FLIP  measurements  were 
computed  according  to: 

%  Difference  =  [  (Xkm  —  Xflip)  /  (Xkm  +  Xflip)  /  2  ]  *  100,  (1) 


1 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
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where  X  is  temperature,  salinity,  ag(Z),  apg(^-)?  cpg(A,),  bbp(^),  y,  or  np  as  a  function  of  depth  and  as  a 
function  of  at  (apg(k)  and  cpg(/_)  are  the  particulate  plus  dissolved  absorption  and  attenuation 
coefficients,  respectively). 

To  first-order,  the  relationship  presented  by  Davies-Colley  (1988)  and  further  examined  by  Zaneveld 
and  Pegau  (2003)  was  used  to  estimate  the  horizontal  visibility  of  a  black  target,  y: 

y  =  4.8/a,  (2) 

where  a  is  the  photopic  beam  attenuation  coefficient,  which  is  a  function  of  the  spectral  background 
radiance  and  the  distance  between  the  target  and  the  observer.  It  was  approximated  by  (Zaneveld  and 
Pegau,  2003): 

a  =  cpg(532)  *  0.9  +  0.081.  (3) 

In  this  regard,  only  one  of  the  IOPs,  the  beam  attenuation  coefficient  at  532  mn,  was  necessary  for 
predictions  of  horizontal  visibility  of  a  black  target.  Equations  (2)  and  (3)  were  applied  to  IOP  data 
collected  during  low  wind,  stratified  conditions  (<  5  m  s'1)  and  high  wind,  mixed  conditions  (>  5  m  s'1) 
and  used  to  examine  the  differences  in  predicted  visibility. 

Objective  (2) 

The  MTF,  which  is  the  magnitude  of  the  Fourier  transform  of  the  PSF,  was  computed  for  two  different 
water  bodies  using  measurements  of  optical  properties  and  analytical  formulations  described  by  Wells 
(1973)  and  Hou  et  al.  (2007).  The  MTF  of  an  aquatic  medium,  H(\j/,  z),  can  be  represented  by: 

H(i|/,z)  =  e'DWZ,  (4) 

where  Z  is  the  range,  D(\|/)  is  the  decay  transfer  function  (DTF),  and  \\i  is  spatial  frequency.  By 
application  of  small-angle  scattering  approximations  (Wells,  1973)  the  collimated  DTF  can  be 
approximated  as: 

Dc(¥)  =  e-S(¥),  (5) 

where  c  is  the  beam  attenuation  coefficient  and  S(\| /)  is  the  light  scattered  back  into  the  acceptance 
cone: 

0max 

Sty)  =  2n  JQ  J0  (2710^)  P(0)  0  d0,  (6) 

where  P(0)  is  the  volume  scattering  function  (VSF),  0  is  the  scattering  angle,  and  J0  is  the  zero-order 
Bessel  function.  The  desired  focused  DTF  can  then  be  obtained  by  transfonnation  in  terms  of  spatial 
frequency  (2ti0\|/)  (see  Wells,  1973  for  more  details): 

D(V)=  l0V(vt)dt.  (7) 

In  situations  where  the  VSF  may  not  be  quantified,  the  following  assumptions  about  the  VSF  have 
been  made  (e.g.,  Wells,  1973): 

P(0)  =  (b  0O)  /  [2ti  (  0O2  +  02) 3/2],  (8a) 

where  0O  is  mean  square  angle  and  b  is  the  scattering  coefficient: 

b  =  2ti  Jq  p(0)  sin  0  d0;  (8b) 

therefore  a  closed-form  solution  for  the  MTF  is  possible  without  direct  measurements  of  the  VSF: 

D(v)  =  c  -  [b  ( 1  -  e'271 00  v)]  /  (2ti  0O  V)-  (9) 

The  optical  data  necessary  for  MTF  computations  were  collected  from  two  field  sites  as  part  of  the 
ONR  sponsored  RaDyO  program:  (1)  Scripps  Institution  of  Oceanography  (SIO)  Pier,  a  shallow-water, 
eutrophic  environment  and  (2)  the  Santa  Barbara  Channel  (SBC),  a  deeper,  mesotrophic  environment 
(Chang  et  al.,  2010).  Profiler  measurements  were  used  to  compute  optical  products:  backscattering 
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ratio,  bb  (k)  =  bhp(A.)/bp(A);  the  slope  of  the  cp(/.)  spectrum,  y,  an  indicator  of  particle  size  distribution 
(Boss  et  ah,  2001);  the  real  part  of  the  index  of  refraction  of  particles,  np,  from  y  and  bh  (k) 
(Twardowski  et  ah,  2001),  useful  for  obtaining  information  about  the  relative  density  of  particles;  the 
single-scattering  albedo,  co0  =  b  /  c;  and  chlorophyll  concentration  (Chi)  from  chlorophyll  fluorescence 
and  spectral  absorption  using  the  methods  of  Sullivan  et  al.  (2005).  The  MASCOT  data  (M. 
Twardowski;  WET  Labs,  Inc.)  were  used  to  compute  a  “bubble  index”  (BI),  where  BI  =  P(70°)  / 
P(120°).  This  BI  is  based  on  theoretically  computed  VSFs  of  relatively  large  bubbles  (>  10  pm),  which 
show  pronounced  enhancement  of  the  VSF  between  60°  and  80°  (Czerski  et  al.,  2011;  Twardowski  et 
al.,  2012).  Hence,  increases  in  the  BI  can  be  indicative  of  a  relative  increase  in  the  concentration  of 
larger  bubbles. 

Relationships  between  physical,  hydrographic,  and  optical  properties  and  the  imaging  perfonnance 
parameter,  MTF,  were  detennined  for  the  SBC  experiment  by  use  of  wavelet  analysis.  Continuous 
wavelet  transforms  (CWTs),  cross  wavelet  transforms  (XWTs),  and  wavelet  coherence  (WTC)  are 
useful  tools  for  analyzing  localized,  intennittent  oscillations  in  data  series,  as  well  as  the  relationship 
and  linkages  between  two  data  series  and  the  strength  of  the  relationship(s)  (see  Grinsted  et  al.,  2004). 
Fourier  expansion  (e.g.,  fast  Fourier  transform;  FFT)  has  traditionally  been  the  technique  employed  for 
detennining  the  frequencies  present  in  a  signal.  However,  FFTs  have  only  frequency  resolution  and  no 
time  resolution.  In  other  words,  although  FFTs  can  reveal  all  frequencies  present  in  a  signal,  they 
cannot  determine  when  these  frequencies  are  present  in  time.  Similar  to  this,  coherence  analysis  can  be 
used  to  find  common  periodicities  between  two  different  time  series;  however,  detennined 
commonality  is  not  localized  in  time  or  space.  Wavelet  analysis  can  be  utilized  to  expand  data  series 
into  time-frequency  space  and  thus  identify  localized  intennittent  periodicities  at  multiple  frequency 
and  temporal  resolutions. 

Objective  (3) 

The  MTF  of  seawater  was  computed  using  profiled  measurements  of  IOPs  conducted  between  the 
near-surface  and  30  m  water  depth  in  SBC  and  Equations  (4)  through  (9)  (see  Chang  et  al.,  2010  and 
Chang  and  Twardowski,  2011).  The  effects  of  stratification  and  scattering  layers  on  the  MTF  were 
examined  for  two  different  scenarios  (Figure  1): 

A.  A  virtual  imager  located  at  the  surface  and  virtual  targets  at  ranges  between  10  cm  and  30  m 
below  the  source. 

B.  Virtual  imagers  located  between  the  surface  and  30  m  and  the  targets  located  10  cm  below  the 
imagers. 
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Figure  1.  Schematic  diagram  of  two  scenarios:  A  (left).  Virtual  imager  (arrow)  and  virtual  targets 
(stars)  separated  by  scattering  layers,  z,„  over  a  total  distance,  Z,  and  B  (right).  Virtual  imagers 

located  a  distance,  z„,  above  virtual  targets. 


For  scenario  A  (Figure  1,  left),  the  collimated  decay  transfer  function  (DTF),  DC(VF),  and  MTF,  H('F), 
at  range,  Z,  and  spatial  frequency,  \|/,  were  computed  following: 

Dc(y)s  =  c  -  z)  dz 

where  the  s  denotes  the  scattering  layer,  c  is  the  beam  attenuation  coefficient  of  layer  s,  S  is  the  light 
scattered  back  to  the  acceptance  cone  of  layer  s,  and  an  integral  taken  between  the  imager  and 
distances  z/,  and  Z2  over  a  total  range  of  Z  to  the  target. 

0max 

S(y)  =  2n  JQ  J0  (27101)/)  P(0)  0  d0 

where  J3(9)  is  the  volume  scattering  function  (VSF),  9  is  the  scattering  angle,  and  J„  is  the  zero-order 
Bessel  function.  The  near-forward  VSF  was  measured  using  a  LISST-100X  and  measurements  of  the 
beam  attenuation  coefficient  were  provided  by  an  ac-s.  The  focused  DTF  was  obtained  by 
transfonnation  in  terms  of  spatial  frequency  (27t0\|/)  prior  to  computation  of  the  MTF  (Wells,  1973): 

D(v|/)=  10'  DC(v|/ 1)  dt 

H(\)/,  z)s  =  e"[Di(v)zi  +  D2(v)z2  + 

Vertical  profiles  of  the  MTF  at  a  spatial  frequency  equal  to  1.0  rad'1  were  constructed  using  the 
equation  for  H(\)/,  z)s  above  by  summing  the  DTFs  between  the  surface  and  each  z„=  10  cm  layer  to  a 
maximum  range  of  Z  =  30  m. 

The  DTF  and  MTF  for  scenario  B  (Figure  1,  right)  were  calculated  following: 

D(V)  =  c  -  S(\|/) 

H(\|/,  z)  =  e'D(v)Z 

where  z  is  the  range  of  each  layer  (e.g.,  zn+2  to  Z),  with  similar  transformation  to  the  focused  DTF.  In 
case  B,  the  vertical  profile  of  the  MTF  for  \|/  equal  to  1.0  rad'1  was  constructed  with  multiple  virtual 
imagers  and  targets,  all  separated  by  the  same  distance  (10  cm)  (Figure  1). 

Scattering  layer  effects  on  the  MTF  for  the  two  different  scenarios  were  quantified  by: 

•  Applying  a  power-law  fit,  y,  to  the  MTF  at  a  given  spatial  frequency  as  a  function  of  depth. 

•  Determining  the  extinction  depth,  k,  of  MTF,  where  the  extinction  depth  was  defined  as  the 
depth  at  which  the  MTF  decreased  to  0.01  cycles/rad. 
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WORK  COMPLETED 


All  work  described  above  was  completed.  Results  are  summarized  below. 

RESULTS 
Objective  (1) 

A  diurnal  wind  pattern  was  observed  during  the  two-week  SBC  experiment.  Winds  were  generally 
calm  (<  4  m/s)  in  the  mornings,  and  increased  to  greater  than  6  m/s,  oftentimes  reaching  10  m/s  by 
1600  PDT.  Starting  on  15  September,  winds  greater  than  about  5  m/s  were  sustained  over  the  course  of 
two  days.  The  persistent  winds  resulted  in  increased  upper  water  column  mixing,  as  evidenced  by  the 
increase  in  30  m  temperature  and  salinity  and  the  deepening  of  the  mixed  layer  depth. 

The  majority  of  optical  variability  between  the  FLIP  and  the  KM  was  due  primarily  to  platform  effects, 
specifically  the  break-down  of  stratification  from  mixing  by  the  hull  of  R/P  FLIP  when  winds  were 
less  than  5  m/s  (Figure  2;  Chang  et  ah,  2010).  Differing  vertical  patterns  in  hydrographic  properties 
and  dissimilar  vertical  distributions  of  water  constituents  (and  therefore  optical  properties)  observed 
between  the  KM  and  FLIP  were  the  consequence  of  FLIP-induced  mixing.  The  FLIP-effect  was  not  a 
factor  during  higher  wind  speeds,  as  the  upper  water  column  was  well-mixed  at  both  platform 
locations. 
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Figure  2.  (A)  Time  series  of  wind  speed.  Crosses  and  triangles  indicate  periods  of  data  collection. 
(B-I)  Profiles  of percent  differences  between  properties  measured  from  the  KM  versus  from  the 
FLIP  as  a  function  of  at The  units  for  ag(412),  apg(412),  cpg(650),  and  bhP(657)  are  m'1.  Crosses 
indicate  periods  of  data  collection  with  wind  speeds  >  5  m/s  and  triangles  represent  wind  speeds  <  5 

m/s  (except  the  last  profile). 

The  observed  variable  vertical  distributions  in  optical  properties  greatly  affected  modeled  underwater 
radiance  distribution  and  predictions  of  horizontal  visibility  (Figure  3).  Hydrolight -modeled 
radiometric  properties  and  AOPs  differed  by  about  50%  between  the  KM  and  FLIP  sites  during 
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stratified,  low  wind  conditions,  as  compared  to  less  than  10%  during  well-mixed  periods.  The 
attenuation-based  model  presented  by  Zaneveld  and  Pegau  (2003)  and  a  backward  Monte -Carlo  model 
were  used  to  predict  the  horizontal  visibility  of  a  black  target  for  high  wind  and  low  wind  conditions. 
Results  from  both  models  were  consistent.  The  results  show  that  highly  variable  IOPs  observed 
between  the  FLIP  and  KM  during  low  winds  can  affect  predictions  of  visibility  by  up  to  57%. 
Differences  were  less  than  15%  between  KM  and  FLIP  predicted  visibility  range  during  well-mixed 
periods. 


Figure  3.  Left :  Vertical  profiles  of  radiometric  and  AOPs  modeled  using  Hydrolight  during  (A-E) 
high  wind  conditions  and  (F-J)  low  wind  conditions.  Pluses  represent  KM  data  and  y  circles  denote 
FLIP  data.  Units  for  Ed  and  Eu  are  W  m~2  mu' ,  for  Lu  are  W  in2  nm'1  sr1,  and  for  r,.s  is  sr1  (R  is 
dimensionless).  Right:  Comparisons  between  computed  horizontal  visibility  of  a  black  target,  y, 
using  data  collected  from  the  KM  versus  from  the  FLIP.  Crosses  and  triangles  denote  data  collected 
during  high  wind  and  low  wind  conditions,  respectively.  Dark  and  light  symbols  signify  results  for 
targets  with  a  diameter  of  0.2  and  0.3  m  (radius  of  0.1  and  0.15  m),  respectively. 


Objective  (2) 

Variability  of  near-surface  MTF  at  both  field  sites  was  strongly  influenced  by  physical  forcing  - 
meteorological  processes  in  the  SBC  and  rip  currents  at  SIO  Pier.  The  2-week  SBC  field  experiment 
can  be  partitioned  into  three  distinct  periods:  (1)  diurnal  winds  (11  -  15  September),  (2)  persistent 
strong  winds  and  upper  water  column  mixing  (15  -  18  September),  and  (3)  an  advection  event 
characterized  by  relatively  high  Chi  and  elevated  optical  quantities  (19-21  September)  (Figure  4). 
During  diurnal  winds  and  the  sustained  strong  wind  event  in  the  SBC,  higher  wind  speeds  were 
associated  with  particle  characteristics  indicative  of  larger,  non-pigmented  particles  with  higher 
relative  indices  of  refraction  and  decreases  in  the  MTF.  The  presence  of  larger  particles  with  higher 
index  of  refraction  observed  during  the  high  winds  of  the  SBC  RaDyO  experiment  could  have  been  the 
result  of  wind-induced  mixing  of  the  upper  water  column  waters  and  surface  bubble  injection  or  an 
increase  in  non-chlorophyll  containing  particles  (perhaps  minerogenic  particles  or  senescent  diatoms) 
or  a  combination  of  both  of  these  factors.  Anderson  et  al.  (2008)  and  references  therein  have  reported 
on  phytoplankton  succession  in  the  SBC,  with  diatom  dominance  during  upwelling  periods  (spring  and 
summer),  followed  by  continued  surface  depletion  of  nutrients  and  stratification,  and  likely  diatom 
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sinking  in  the  autumn.  It  is  possible  that  these  senescent  diatoms  were  retained  in  near  surface  waters 
during  strong  mixing  associated  with  high  winds. 
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Figure  4.  Time  series  of  (a)  wind  stress  (blue)  and  wind  speed  (red),  (b)  cpg(532)  (brown  dots)  and 
Chi  (green  triangles)  at  2  m,  (c)  near-surface  np  (blue  dots)  and  /(red  triangles),  near-surface  (d) 
volumetric  bubble  concentration  derived  from  BI,  and  (e)  MTF  for  W=  3  rad1  (black  dots)  and  *F= 
50  rad'1  (gray  triangles)  calculated  from  2  m  optical  data  collected  in  the  SBC. 

Surface  bubble  injection  is  often  associated  with  strong  sustained  winds,  which  in  itself  would  result  in 
substantive  changes  in  particle  composition  characteristics  in  surface  waters  (Zhang  et  ah,  2011). 
Information  about  bubble  concentration  was  obtained  from  the  BI,  which  was  derived  from  the  ratio  of 
MASCOT -measured  VSFs  at  70°  and  120°.  The  BI  can  be  used  as  a  proxy  for  the  concentration  of  the 
large  bubble  subpopulation,  following  Twardowski  et  al.  (2012).  The  derived  concentration  of 
relatively  large  bubbles  (>  10  pm  in  diameter)  increased  by  nearly  an  order  of  magnitude  on  15 
September  2008  and  by  a  factor  of  three  on  18  September,  coinciding  with  periods  of  sustained 
elevated  wind  speeds  (Figure  4d).  The  increase  in  the  concentration  of  the  large  bubble  subpopulation 
was  accompanied  by  a  decrease  in  the  slope  of  the  particulate  attenuation  spectra,  indicative  of  an 
increase  in  particle  size  (Figure  4c). 

Although  not  statistically  significant  in  2-m  time  series  data,  the  observed  advection/bloom  event 
starting  on  19  September  2008  resulted  in  significant  coherence  between  the  MTF  and  the  IOPs  and 
IOP  products  indicative  of  particle  concentration  and  composition  at  about  20-m  water  depth.  This 
event  was  marked  by  a  decrease  in  the  MTF,  increases  in  the  magnitude  of  all  IOPs,  Chi,  and  np  and  a 
decrease  in  y,  indicative  of  an  increase  in  the  concentration  of  larger  biogenic  particles  with  relatively 
higher  refractive  index  (Figure  4).  The  concentration  of  the  large  bubble  subpopulation  during  this 
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time  period  was  relatively  low,  suggesting  lesser  effects  due  to  surface  bubble  injection  and  more  to  a 
phytoplankton  bloom  that  was  associated  with  a  shift  in  the  SBC  circulation  (Figure  4). 

Highly  reflective  minerogenic  particles  with  relatively  high  bulk  particle  index  of  refraction  were 
associated  with  SIO  Pier  rip  currents,  which  resulted  in  a  large  increase  in  the  MTF.  These  particles, 
although  not  explicitly  sampled  or  analyzed,  appeared  like  mica  flakes,  which  although  rare  in  beach 
sand,  have  been  observed  along  beaches  in  San  Diego,  CA.  This  increase  in  the  MTF  was  likely  due  to 
enhanced  reflectivity  coupled  with  relatively  high  visibility;  beam  attenuation  was  not  greatly  elevated 
during  the  rip  current  event  and  Chi  was  at  a  minimum  (Figure  5). 
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Figure  5.  SIO  Pier  time  series  of2-m  (a)  cpg(532)  (brown  dots)  and  Chi  (green  triangles)  and  (b) 
bhP(532)  (blue  dots)  and  np  (red  triangles);  and  (c)  MTF  (tF=  3  rad'1  in  black  dots  and  XF=  50  rad1 

in  gray  squares). 

These  observations  are  particularly  interesting  as  one  would  expect  the  imaging  performance 
parameter  to  be  influenced  primarily  by  particle  concentration  as  revealed  by  the  magnitude  of  IOPs; 
however  it  was  found  that  particle  size  and  type  (relative  density  information  gained  from  the  real 
index  of  refraction  of  particles)  were  important  sources  of  variability  in  the  MTF. 

The  results  show  that  the  effects  of  particle  composition  on  the  shape  of  the  VSF  in  the  near-forward 
direction  are  important  to  imaging  performance,  which  is  a  unique  observation  facilitated  by 
advancements  in  in  situ  optical  instrumentation.  By  performing  a  sensitivity  analysis  using  the  [!- 
approximation  formulation  for  the  MTF,  the  relative  importance  of  particle  effects  on  the  MTF  was 
determined  (Figure  6).  Results  suggest  that  at  0O  between  0.07  and  0. 1 1  rad,  particle  concentration  and 
composition  affected  MTF  variability  similarly.  At  90  greater  than  0.1 1  rad,  MTF  variability  was 
stronger  with  variable  IOPs,  i.e.  particle  concentration  effects  dominated.  Strikingly,  when  90  is 
between  0.03  and  0.07  rad  (i.e.  particle  size  distributions  trending  toward  larger  sizes),  the  MTF  is 
much  higher  for  frequencies  less  than  30  rad"1  (Figure  6).  Here,  particle  composition  effects  dominated 
MTF  variability,  in  a  relative  sense.  For  the  specific  ranges  of  variability  in  optical  properties  observed 
during  the  study,  these  modeling  results  are  consistent  with  the  strong  coherence  observed  between 
particle  composition  parameters  and  imaging  performance. 
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Figure  6.  Particle  concentration  versus  composition  effects  on  MTF  variability,  (a)  MTF  computed 
by  varying  the  beam  attenuation  coefficient,  c  (as  indicated),  while  keeping  the  mean  square  angle, 
90,  constant  at  0.07  rad.  (b)  MTF  computed  by  varying  G„  (as  indicated)  while  keeping  c  constant  at 
0.6  m~.  The  range  was  set  at  2  m  and  a>„  was  0.9  in  both  situations. 

The  small-angle  approximation  (Eq.  8;  e.g.,  Wells,  1973)  assumes  strong  forward  scattering,  which 
was  not  the  case  during  the  SIO  Pier  rip  currents.  To  explore  this  further,  MTFs  computed  from  the 
DTF  using  measured  P  (Eq.  7)  were  evaluated  versus  those  calculated  using  approximated  P  (Eq.  9;  90 
=  0.07  rad)  for  the  SIO  Pier  data  set  (Figure  7).  The  p -approximation  method  produced  very  similar 
MTFs  to  those  computed  using  measured  p,  particularly  at  high  spatial  frequencies  (>  10  rad’1).  MTFs 
at  spatial  frequencies  less  than  10  rad"1  were  not  comparable  during  periods  when  the  backscattering 
coefficient  exceeded  0.012  m"1.  Results  from  the  two  methods  deviated  greatly  from  each  other  at  all 
spatial  frequencies  during  the  rip  current  event  (Figure  7).  Here,  it  is  shown  that  in  highly 
backscattering  environments  such  as  those  with  substantial  minerogenic  particles,  the  shape  of  the  VSF 
cannot  be  assumed  when  computing  imaging  performance  parameters  such  as  the  MTF. 
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Figure  7.  SIO  Pier  time  series  of  MTF  at  x¥=  (a)  3  rad'1  and  (b)  50  rad'1,  with  MTFs  computed 
using  the  /3-approximation  (e.g.,  Wells,  1973)  shown  with  gray  triangles. 
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Objective  (3) 

The  upper  water  column  MTF  at  small  spatial  frequencies  in  the  presence  of  scattering  layers  and 
stratification  can  be  orders  of  magnitude  less  than  the  MTF  in  mixed  conditions,  as  illustrated  in  Figure 
8.  The  distance  between  the  imager  and  target  contributed  greatly  to  the  magnitude  of  the  MTF  when 
imager  and  target  were  separated  by  a  scattering  layer. 


MTFpP  =  1)  MTF(T  =  1) 


Figure  8.  Vertical  profiles  of  the  MTF(X¥=  1  rad1)  shown  to  a  range  of  Z  =  15  m  for  A:  Scenario  A 
and  B.  Scenario  B.  The  solid  dots  indicate  the  extinction  depth,  tc  Profile  color-coding  denotes 
specific  periods  of  the  SBC  RaDyO  experiment:  blues  =  strong  stratification,  yellow-green  =  upper 
water  column  mixing,  and  reds  =  strong  scattering  layers  from  localized  upwelling  (see  Chang  et  al., 
2010  and  Chang  and  Twardowski,  2011  for  more  details). 

The  extinction  depth  determined  for  the  MTF(\j/  =  1  rad'1)  with  an  imager  at  the  surface  (scenario  A) 
did  not  vary  greatly  from  about  4.5  m  for  all  15  profiles;  the  minimum  value  of  k  was  2.9  m  in  the 
presence  of  strong  scattering  layers  and  maximum  k  was  4.8  m  during  mixed  conditions  (Figure  9). 
When  virtual  imagers  and  targets  were  placed  10  cm  apart  throughout  the  water  column,  k  varied 
between  4.5  m  (strong  scattering  layers)  and  9.8  m  (mixed  conditions). 

The  differences  in  y  between  scenarios  A  and  B  were  larger  during  stratified  conditions  (average 
difference  of  1.5)  as  compared  to  during  upper  water  column  mixed  conditions  (average  difference  of 
1.25).  The  smallest  difference  in  y  between  scenarios  A  and  B  (0.985)  occurred  during  periods  with 
strong  scattering  layers  (Figure  9). 
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Figure  9.  Upper:  MTF('F=  1  rad1)  extinction  depth,  Middle:  Slope  of  the  power-law  fit  to  the 
vertical  profile  of  MTF(XF=  1  rad1),  and  Lower:  The  mixed  layer  depth  computed  as  the  depth  at 
which  the  cpg(660)  change  from  near-surface  cpg(660)  is  greater  than  0.05  in1.  Yellow  =  scenario  A 
and  Purple  =  scenario  B.  The  profile  numbers  can  be  characterized  by:  1-7  =  strong  stratification, 

8-10  =  mixed,  and  11-15  =  strong  scattering  layers. 


These  results  imply  that  in  the  presence  of  near-surface  optical  scattering  layers,  an  underwater  image 
taken  from  just  below  the  surface  to  a  depth  of  a  mere  0.5  m  can  be  significantly  blurrier  than  if  the 
imager  was  located  tens  of  cm  below  the  surface.  Additionally,  analysis  or  reconstruction  of  an  image 
taken  by  a  near-surface  detector  of  a  target  located  at  just  a  couple  of  meters  below  the  surface  may  not 
be  possible  in  the  presence  of  strong  scattering  layers,  even  if  the  water  appears  relatively  clear  from 
the  surface. 

IMPACT/APPLICATIONS 

Objective  (1) 

The  results  have  the  potential  to  greatly  impact  the  modeling  of  underwater  radiance  distribution  for 
the  RaDyO  project.  The  hull  of  FLIP  itself  has  the  capability  to  mix  the  upper  water  column  and  affect 
SBL  processes.  Physical  measurements  conducted  by  RaDyO  Pis  off  the  FLIP  could  have  been 
compromised.  Additionally,  IOPs  were  not  collected  off  the  FLIP  during  the  Hawaii  experiment;  hence 
the  model(s)  would  have  been  constructed  and/or  validated  with  physical,  radiometric,  and  apparent 
optical  property  (AOP)  measurements  collected  off  the  FLIP  and  IOP  measurements  collected  off  the 
KM,  which  was  several  km  away. 

Objective  (2) 

Optically-derived  particle  characteristics  such  as  bulk  particle  index  of  refraction  and  particle  size 
distribution  were  shown  to  be  significantly  related  to  the  variability  of  the  MTF.  This,  in  addition  to 
the  great  temporal  variability  in  environmental  conditions  observed  at  two  different  field  sites  stresses 
the  importance  of  obtaining  or  predicting  optical  properties  and  derived  particle  composition 
characteristics  (describing  relative  particle  size  and  type)  at  temporal  and  spatial  scales  necessary  to 
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quantify  rapid  changes  in  image  processing  parameters.  Importantly,  estimates  of  water  clarity  or 
turbidity  are  not  the  only  factors  influencing  the  variability  of  the  image  performance  parameters, 
particularly  at  steep  near-forward  scattering  angles,  where  variability  in  the  MTF  was  found  to  be 
dominated  by  particle  composition  effects.  Additionally,  the  standard  approximation  for  the  VSF 
cannot  be  used  to  compute  imaging  performance  in  environments  with  highly  backscattering/reflective 
particles. 

Objective  (3) 

Knowledge  of  the  MTF  is  important  for  the  interpretation  of  images  from  underwater  electro-optical 
imaging  systems.  MTF-based  image  analysis  and  reconstruction  is  influenced  by  optical 
inhomogeneity  (i.e.  scattering  layers)  as  well  as  by  the  vertical  placement  of  the  target  versus  that  of 
the  detector,  or  imager  (Wells,  1973;  Frew  and  Voss,  1997).  Take,  for  example,  the  remark  made  by 
Wells  (1973),  “You  can  lay  a  sheet  of  vellum  on  a  printed  page  and  read  through  it  quite  well,  but  the 
same  vellum  lifted  a  centimeter  off  the  page  completely  obscures  vision.”  This  statement  is  now 
outdated  as  vellum  is  no  longer  a  common  source  for  paper.  Here  is  another  way  to  consider  the 
problem:  one  should  not  feel  embarrassed  while  showering  several  feet  away  from  a  textured  shower 
door;  however  body  parts  could  be  fully  revealed  should  they  make  contact  with  the  door  (Figure  10). 


Figure  10.  These  photographs  of  Peso  the  dog  were  manipulated  in  Photoshop  to  simulate  Peso  well 
behind  a  textured  glass  shower  door  (Left)  and  Peso  ’ s  nose  pressed  against  the  glass  (Right). 

This  same  effect  applies  in  stratified  ocean  environments  where  a  target  and  imager  are  separated  by 
optically  distinct  layers.  These  layers  could  have  important  implications  on  imaging  performance  and 
analysis.  The  ability  to  determine  optimal  placement  of  passive  imagers  relative  to  target  locations 
would  constitute  dramatic  improvements  to  image  clarity  and  restoration  efforts.  Additionally,  image 
reconstruction  would  be  facilitated  with  knowledge  of  the  effects  of  scattering  layers  on  image 
distortion  and  imaging  performance  parameters. 

RELATED  PROJECTS 

At  present,  none. 
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